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Abstract 

A code has been developed to automatically sim- 
plify full chemical mechanisms. The method em- 
ployed is based on the Intrinsic Low Dimensional 
Manifold (ILDM) method of Maas and Pope 1,2 . 

The ILDM method is a dynamical systems ap- 
proach to the simplification of large chemical kinetic 
mechanisms. By identifying low-dimensional at- 
tracting manifolds, the method allows complex full 
mechanisms to be parameterized by just a few vari- 
ables: In effect, generating reduced chemical mech- 
anisms by an automatic procedure. These resulting 
mechanisms however, still retain all the species used 
in the full mechanism. 

Full and skeletal mechanisms for various fuels are 
simplified to a two dimensional manifold, and the 
resulting mechanisms are found to compare well 
with the full mechanisms, and show significant im- 
provement over global one step mechanisms, such 
as those by Westbrook and Dryer 3 . In addition, 
by using an ILDM reaction mechanism in a CFD 
code, a considerable improvement in turn-around 
time can be achieved. 

Introduction 

Detailed chemical mechanisms for hydrocarbon 
fuels typically contain in order of 50 species and 200 
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rate equations. Applying these mechanisms to a 3D 
turbulent combustion problem is beyond the capa- 
bilities of todays computers, and beyond those in 
the foreseeable future. Thus the practice of simpli- 
fying the full mechanisms is used, where the number 
of species and equations is reduced by various as- 
sumptions about the reaction, e.g. Partial equilib- 
rium of species. In the crudest sense, this reduction 
is achieved by throwing out, or combining species 
and rate equations. 

However to perform these simplifications, consid- 
erable human input is required, both in knowledge 
of the detailed mechanism and testing the resulting 
reduced mechanism. In addition, the accuracy of 
the reduced mechanism drops with every simplifi- 
cation, requiring the user to compromise between 
how much accuracy they need against how much 
they can afford. 

This need for considerable human effort to ob- 
tain reduced mechanisms was the motivation for 
the development of the ILDM method 1,2 of sim- 
plifying chemical kinetics. The objectives were to 
provide a completely automated method of simpli- 
fying chemical kinetics, with the only input being 
the full mechanism, and the degree of simplification. 

In this paper, the development of an ILDM code 
will be described. This code will be used to simplify 
several hydrocarbon mechanisms, and the perfor- 
mance of these simplifications will be tested in a 
Perfectly Stirred Reactor (PSR) code. The ILDM 
results will be compared against the full mechanism, 
a 12 species, 10 step reduced mechanism and a one 
step mechanism. In addition, the performance of 
the reduced mechanism in a CFD code will be in- 
vestigated. 
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ILDM Method 

The composition of a chemical system can be 
thought of as a point in a multidimensional space: 
each dimension representing a species concentra- 
tion. Thus the rate equations can be thought of 
as representing how that point moves in the scalar 
space. The ILDM approach is based on analysing 
and simplifying these equations of motion of the 
scalar point. 

How the ILDM method does this is to identify 
attracting manifolds in the scalar space. Because 
of the wide range of time scales in a chemical sys- 
tem, the composition changes very rapidly in some 
directions and moves slowly in others. Thus the 
composition moves rapidly toward the attracting 
manifold, and then moves along that surface to 
the equilibrium point. To simplify the mechanism, 
one simply assumes the composition always lies on 
these manifolds. Thus all one need do is specify 
the dimension of the manifold. For example, a 
zero-dimensional manifold is simply the equilibrium 
point while a one-dimensional manifold is a line, on 
which lies the equilibrium point. 

Because a one-dimensional manifold is just a line 
in the scalar space, it can be parameterized by just 
one variable. This variable can be a scalar, a combi- 
nation of scalars or even a thermodynamic property, 
and is referred to as the reference variable. The only 
requirement is that the reference variable does not 
give a multi-valued function. 

One interesting benefit of this method of simpli- 
fication is that the value of all the species in the full 
mechanism on the manifold are available, as a func- 
tion of the reference variable. Thus the minor, and 
trace species are available without any increase in 
the complexity of the reaction mechanism. This is 
contrasted with conventional reduced mechanisms 
that require additional species and equations to be 
added if more detailed composition information is 
required. In addition, state properties can also be 
stored on the manifold, and so properties like den- 
sity are calculated with all the minor species being 
taken into account , thus providing a more accurate 
value. 

In the code developed at Lewis, the ILDM sim- 
plified mechanism is obtained by a trajectory- 
generated technique, described by Pope and Maas 
4 . In this method, the full mechanism is param- 
eterized by two scalars, a mixture fraction based 
on atom concentrations and a progress variable, 
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based on either species mass fractions, temperature 
or Gibbs function. The resulting species concen- 
trations, rates and properties are stored in look-up 
tables. Also provided is a set of interpolation rou- 
tines that can be called from CFD codes. These 
provide the interface between the ILDM tables and 
the combustion code. 

Reaction Mechanisms 

In this paper, three hydrocarbon fuels are stud- 
ied: methane, propane and Jet A. The mechanisms 
used are the methane mechanism of Miller and 
Bowman 5 , the propane mechanism of Jachimowski 
6 and the Jet A mechanism of Kollrack ' . 

In addition, two reduced mechanisms are used 
for comparison. These are the one step methane 
mechanism of Westbrook and Dryer 3 and the 12 
step methane mechanism of Kundu 8 . 

Test Case 

The test case employed for the reaction kinet- 
ics is the Perfectly Stirred Reactor (PSR). In this 
example, a mixture of fuel and air flows into a re- 
actor vessel. This mixture mixes instantly with the 
composition in the container, reacts and the result- 
ing products exit the container. Residence time is 
defined as the density of the exit composition mul- 
tiplied by the volume of the reactor, divided by the 
mass flow rate of the system. 

Initial conditions are chosen to be those similar 
to gas turbine operation, with fuel and air inlet 
temperatures at 350K and 750K respectively, and a 
pressure of 10 atmospheres. 

Results 

Results will consist of plots of residence time vs. 
property for the several different reaction schemes 
and for different inlet mixture fractions. Due to the 
quantity of results available, only a selection will be 
presented. In all these cases, the progress variable 
( Yp) chosen is the sum of the mass fractions of CO 2 
and H20. 

The first area addressed is the table resolution, 
i.e. How many grid points are required to accu- 
rately represent the reaction scheme. A stoichio- 
metric mixture of methane and air was chosen as 
the representive composition for this test. 
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In Fig. 1, the results for refinement of the table in 
the Yp direction are shown against the full mecha- 
nism. It can be seen that having 32 points provides 
very good agreement with the full mechanism for 
all species. 

In Fig. 2, the effect of table resolution in the 
mixture fraction (xi) direction is investigated. It 
can be seen that a grid spacing of 0.0025 provides a 
very good agreement with the full mechanism. This 
corresponds to 24 grid nodes between the flamibility 
limits. Considering these results, a grid resolution 
of 32x32 is used for the remaining cases. 

In Fig. 3 and 4, the performance of the ILDM 
table is shown compared to the full mechanism 5 , 
a 12 species, 10 step mechanism 8 and a one step 
mechanism 3 . It can be seen that the ILDM scheme 
compares very well with the full mechanism and the 
10 step scheme in all plots. Especially noteworthy 
is the very good agreement with the density, as this 
variable is the one that affects the flow calculations. 
The poor agreement by the one-step mechanism can 
be attributed to the lack of minor species in that 
mechanism. It should be noted however that the 
results for the mass fraction of fuel show a very 
poor agreement. This lack of agreement is present 
in all the cases tested, and will be discussed in the 
next section. 

The same fuel was then used, but for lean and 
rich mixtures. The results of these are shown in 
Fig. 5. Both give good agreement with the full 
mechanism. 

In Fig. 6, the results for a stoichiometric mixture 
of propane and air are shown. Again the results are 
in good agreement with the full mechanism. The 
only difference observed is the kink in the Yh 20 
result for small residence times. However the very 
good agreement of the density result shown that 
this would not have a noticeable effect on the overall 
performance of the table. 

However the results for JetA, shown in Fig. 7 
show a divergence from the full mechanism results 
at short residence times. This can be probably be 
attributed to the full mechanism used in this test. 
The mechanism used 7 is the smallest of the three 
mechanisms, despite being the most complex fuel. 
The problem may lie in the omission of some impor- 
tant intermediate species. In addition, this mecha- 
nism is the oldest of all three. Comparing the ele- 
mentary rate constants with those of methane and 
propane reveal many differences. It is quite pos- 
sible that the rate constants of some intermediate 
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species may be in error, thus affecting the high-load 
results. 

The next part of the mechanism investigated was 
the speed of the look-up table compared to using 
an equivalent reduced mechanism. To perform this 
test, the 12 species methane mechanism used in a 
combustion code was analyzed. For a 500,000 node 
problem, the CFD run-time was about 14.0 seconds 
per iteration, using 84 nodes on an SP2 machine. 
Of that time, 7.0 seconds was spent in performing 
the reaction and calculating the composition prop- 
erties. 

In this code, the equivalent look-up table require- 
ments would require 0.17 sec. In addition, because 
the transport equations for only mixture fraction 
and progress variable would need to be calculated, 
rather than 12 species, it is estimated that the time 
saved using the lookup tables would be about 10.7 
seconds per iteration. 

The space required for these look-up tables is es- 
timated as 4.3KBytes per property stored. This 
gives a total table size of about 60KBytes. 

Discussion 

A feature of the results that needs to be discussed 
is the poor prediction of fuel concentrations in all 
three cases. 

For the fuel prediction, the lack of agreement 
between the full mechanism results and the re- 
duced mechanisms is a function of the simplifica- 
tion procedure. Put simply, the break-up of the 
fuel molecule is one of the processes that the ILDM 
scheme deems to occur very rapidly, and so the 
PSR results show essentially zero concentration un- 
til near extinction. However the performance of 
the model, specifically in the prediction of density, 
shows that the error in predicting the fuel concen- 
trations does not compromise the overall accuracy 
of the model. 

Another area to be addressed is the implementa- 
tion of the ILDM mechanism in a CFD code. 

The current methodology used is to store the 
mechanism in a simple look-up table as a function of 
the mixture fraction and the progress variable. The 
look-up table consists of a uniform grid, bounded in 
the mixture fraction direction by the lean and rich 
flamability limits. In the progress variable direction 
it is bounded by the pure mixing and fully burnt 
limits. When used in a CFD code, it is found that 
the best way to implement the scheme is by storing 
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the reaction increments for a set of different time 
steps. In addition the properties such as density, 
laminar viscosity and thermal diffusivity are also 
stored in the table. Thus the interpolation routine 
is given an initial composition and time step, and 
it returns the updated composition, along with the 
new property values. 

It should also be noted that the current imple- 
mentation of the ILDM scheme contains some lim- 
itations to its applicability. The assumptions made 
for the current model are constant pressure, unity 
Lewis number and simple inlet conditions, i.e. 
All fuel inlet temperatures and compositions must 
be the same, and all air inlet temperatures must 
have their temperatures and compositions identi- 
cal. However these are just restrictions of the cur- 
rent implementation, not of the ILDM scheme in 
general. By including pressure as a tabulation vari- 
able, the uniform pressure restriction can be lifted. 
By including other variables such as enthalpy and 
composition the other cases can also be solved for 
by this method. 

Conclusion 

In this paper it has been shown that a 1-step 
ILDM scheme gives good agreement with full mech- 
anism results in the PSR test for a variety of hy- 
drocarbon fuels. The accuracy of the method has 
also been shown to be similar similar to that of a 
12 species, 10 step mechanism. In addition due to 
the tabular nature of the storage, the ILDM scheme 
can provide significant savings in CPU time when 
coupled with a CFD code. While the current im- 
plementation does restrict the class of problem that 
can be treated by this method, it should be noted 
that relatively simple changes to the storage meth- 
ods can overcome these. 
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Figure 1: Test of table resolution in Y P direction for stoichiometric methane-air mixture. PSR results 
shown as composition quantities plotted against residence time. O " Full mechanism; □ - 64 nodes; A - 
32 nodes; <1-16 nodes; V " 8 nodes. 
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Figure 2: Test of table resolution in xi direction for stoichiometric methane-air mixture. PSR results shown 
as composition quantities plotted against residence time. 0 - Full mechanism; □ - 0.0025 spacing; A - 
0.005 spacing; <1 - 0.01 spacing; y - 0.02 spacing. 
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Figure 3: Comparison of different mechanisms for stoichiometric methane-air mixture. PSR results shown 
as composition quantites plotted against residence time. Q “ Full mechanism; □ - ILDM table; A - 
12-species, 10-step mechanism; <1 - 0.01 Westbrook and Dryer 1-step mechanism. 
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Figure 4: Comparison of different mechanisms for stoichiometric methane-air mixture. PSR results shown 
as composition quantites plotted against residence time. 0 - Full mechanism; □ - ILDM table; A - 
12-species, 10-step mechanism; <3 - 0.01 Westbrook and Dryer 1-step mechanism. 
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Figure 5: Performance of the ILDM mechanism for lean and rich mixtures of methane-air. PSR results 
shown as composition quantites plotted against residence time. Q - Full mechanism; □ - ILDM table. 
Lean mixture with £=0.04 shown in upper row. Rich mixture with £=0.07 in lower row. 
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Figure 6: Performance of the ILDM mechanism for stoichiometric mixture of propane-air. PSR results 
shown as composition quantites plotted against residence time. 0 - Full mechanism; □ - ILDM table. 
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Figure 7: Performance of the ILDM mechanism for stoichiometric mixture of Jet A-air. PSR results shown 
as composition quantites plotted against residence time. 0 - Full mechanism; □ - ILDM table. 


NASA CR-204138 


11 







REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson 
Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188), Washington, DC 20503 


1. AGENCY USE ONLY ( Leave blank) 


4. TITLE AND SUBTITLE 


2. REPORT DATE 

August 1997 


Automated Simplification of Full Chemical Mechanisms 



3. REPORT TYPE AND DATES COVERED 

Contractor Report 


5. FUNDING NUMBERS 


WU-523-26-33-00 

NCC3-534 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

Institute for Computational Mechanics in Propulsion 
22800 Cedar Point Road 
Cleveland, Ohio 44142 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


E-10855 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESSEES) 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 


11. SUPPLEMENTARY NOTES 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 

NASA CR-204138 
I COMP-97 -07; CMOTT-97-02 
AIAA-97-31 15 


ICOMP Program Director, Louis A. Povinelli, organization code 5800, (216) 433-5818. 


12a. DISTRIBUTION/AVAILABILITY STATEMENT 


12b. DISTRIBUTION CODE 


Unclassified - Unlimited 
Subject Category 34 

This publication is available from the NASA Center for AeroSpace Information, (301) 621-0390. 


13. ABSTRACT ( Maximum 200 words) 

A code has been developed to automatically simplify full chemical mechanisms. The method employed is based on the 
Intrinsic Low Dimensional Manifold (ELDM) method of Maas and Pope. The ILDM method is a dynamical systems 
approach to the simplification of large chemical kinetic mechanisms. By identifying low-dimensional attracting mani- 
folds, the method allows complex full mechanisms to be parameterized by just a few variables: In effect, generating 
reduced chemical mechanisms by an automatic procedure. These resulting mechanisms however, still retain all the 
species used in the full mechanism. Full and skeletal mechanisms for various fuels are simplified to a two dimensional 
manifold, and the resulting mechanisms are found to compare well with the full mechanisms, and show significant 
improvement over global one step mechanisms, such as those by Westbrook and Dryer. In addition, by using an ILDM 
reaction mechanism in a CFD code, a considerable improvement in turn-around time can be achieved. 


14. SUBJECT TERMS 



Combustion; Kinetics; Reduced mechanisms; Manifold methods 

17. SECURITY CLASSIFICATION 

18. SECURITY CLASSIFICATION 

19. SECURITY CLASSIFICATION 

OF REPORT 

OF THIS PAGE 

OF ABSTRACT 

Unclassified 

Unclassified 

Unclassified 


15. NUMBER OF PAGES 

13 


16. PRICE CODE 

A03 



NSN 7540-01-280-5500 


Standard Form 298 (Rev. 2-89) 

Prescribed by ANSI Std. Z39-16 
29B-102 




















